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Abstract 

It is shown that most of the available data on the pp — > ppr] reaction, including the invariant 
mass distributions in the pp pprj reaction recently measured at COSY, can be understood in 
terms of the partial-wave amplitudes involving final pp S and P states and the r] meson s-wave. 
This finding, together with the fact that results within a meson-exchange model are especially 
sensitive to the details of the excitation mechanism of the /S'ii(1535) resonance, demonstrates the 
possibility of investigating the properties of this resonance in A^A^ collisions. The spin correlation 
function Cxx is shown to disentangle the S- and P-wave contributions. It is also argued that spin 
correlations may be used to help constrain the contributions of the amplitudes corresponding to 
the final pp '^Pq and '^^2 states. 

PACS numbers: PACS: 25.10.+S, 13.75.-n, 25.40.-h 
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The primary motivations for studying the production of mesons off nucleons and nuclei 
are to investigate the structure and properties of the nucleon resonances and to learn about 
hadron dynamics at short range. As far as hadron-induced reactions are concerned, and 
specifically nucleon-nucleon {NN) collisions, there is already a wealth of information on 
the production of the lightest meson, the pion. In particular, there now exists a fairly 
accurate and complete set of data, especially for tt° production in the near-threshold energy 
region Q|, which should allow for a partial wave analysis. The 1] meson, which is the 
next lightest non-strange meson in the meson mass spectrum, has also been the subject 
of considerable interest. A peculiar feature of this meson is that it couples strongly to 
the S'ii(1535) nucleon resonance, which offers a unique opportunity for investigating the 
properties of this resonance. Unfortunately, the experimental information on t] production 
in NN collisions 0, S, 0, S, 0, Q] is much less complete than for pion production and is not 
yet sufficient for a model-independent partial-wave analysis. However, the available data 
base has greatly expanded recently thanks to measurements by the TOF and COSY-11 
collaborations at COSY 0, Q] that provided, not only r] and proton angular distributions, 
but also pp and pr] invariant mass distributions for the reaction pp — > ppr]. These new data, 
together with the earlier measurements , open the possibility for investigating this 

reaction in much more detail than could be done previously. 

A general feature of meson production in A^A^ collisions is that the energy dependence 
of the total cross section in the near-threshold region is basically dictated by the available 
phase space plus the A^A^ final state interaction (FSI) in S states. The effect of the strong 
A^A^ FSI also shows up in the corresponding A^A^ invariant mass spectrum as a peak close 
to the threshold value of the invariant mass, m^vAr = 2m jv, where denotes the nucleon 
mass. Surprisingly, the recently measured pp invariant mass distribution in the reaction 
pp ppr] p. Ill at excess energies of Q = 15 and 41 MeV shows, in addition to a peak very 
close to the threshold, a broad bump at higher values of rripp (see Fig. H}. While the peak 
can easily be understood as arising from the strong pp FSI in the ^5*0 state as mentioned 
above, it is not trivial to explain the origin of the bump at higher rupp. 

The purpose of the present work is to analyze this seemingly peculiar feature exhibited 
by the pp invariant mass distribution. Thereby we will show that the available data on 
pp —>■ ppr], including the invariant mass distribution, can essentially be understood in terms 
of S- and P-wave amplitudes. This result suggests that the properties of the 5*11(1535) 
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nucleoli resonance can be studied here in terms of only a few partial-wave amplitudes. We 
start by considering the possible partial-wave states in pp —>■ pprj near the threshold j^. 
First, in this reaction, the t] meson is dominantly produced through the excitation and de- 
excitation of the 5*11(1535) resonance. Therefore we expect that the 1] meson should be 
produced mainly in the s-wave. Of course, this is strictly true only in the rest frame of 
the resonance and not necessarily in the overall cm. frame of the system. However, near 
threshold, this should not make a significant difference. In fact, the observed 77 angular 
distribution in the overall cm. frame P| is practically isotropic. In addition, the very first 
analyzing power measurement in pp pprq by the COSY-11 group yielded rather small 
values. Indeed this observable is basically consistent with zero, given the relatively large 
uncertainties involved, and therefore consistent with pure s-wave contributions. As long 
as we restrict ourselves to an ri meson in the s-wave and final A^A^ state to the S and P 
wave, we have o.lv f^ee pa.t,a.-wave a„,p.,tude. that can eon.ibute to this .action fl: 
■^Po ^Sqs, ^Sq ^Pqs and ^P2S. Among these, we would naively expect that 

the ^Pq —>■ ^Sqs is the only relevant contribution near threshold. However, as mentioned 
above, the contribution of the S'-wave alone is unable to explain the observed bump in the 
pp invariant mass distribution. 

One plausible explanation may be attributed to effects from the riN FSI. Indeed there 
are already strong indications from the total production cross sections that the rjN FSI may 
play an important role in the reaction A^A^ —>■ NNrj near threshold. For both pp and pn 
induced rj productions one has observed that there is an enhancement of the production 
cross sections for very small excess energies that cannot be explained by the A^A^ FSI effects 
alone 0,0]. However, those effects seem to be confined to an excess energy range of up to 
at most 20 MeV from the threshold. Thus, one would expect that the rjN FSI effects should 
have an influence on the invariant mass spectrum measured at the lower energy of Q = 15 
MeV. It is, however, unlikely that such effects should still be so important at Q = 41 MeV. 
A proper inclusion of the riN FSI calls for solving the Faddeev equation in the three-particle 
continuum which is technically very involved. A rough estimate suggests that a rather strong 
riN interaction would be needed to reproduce the data at the higher energy which is 
difficult to be reconciled with other information about the r]N interaction. Therefore, we 
seek an alternative explanation based largely on the observation that the shape of the bump 
can be reproduced by folding p'"^ with the available phase space. Here, p' denotes the relative 
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momentum in the final pp system. This suggests to us, that the bump seen in the experiment 
could be simply caused by the pp P-wave in the final state. Admittedly, the measured final 
proton angular distribution in the overall cm. frame is nearly isotropic 0], which could be 
seen as an evidence against large P-wave and higher partial wave contributions. However, 
as we shall show below, there is no principal contradiction between a nearly isotropic proton 
distribution and a significant P-wave fraction in the invariant mass spectrum. 

Let us now make some general remarks about the structure of the reaction amplitude 
for pp pprj. In what follows, we shall assume that the rj is in an s wave relative to the 
(final) pp system and that the final protons are in a relative S and/or P state. Since the r) 
meson is an isoscalar pseudoscalar particle, it follows immediately that the orbital angular 
momentum of the pp system has to change in the transition from the initial to the final 
state and consequently, due to the Pauli principle, the total spin also has to change. Thus, 
the most general form of the reaction matrix (involving even angular momenta of rj) can be 
written as 

M = (aP5'=o - BPs'=i) -^i^i- (^2) , (1) 

where Ps'=o,i stands for the total spin singlet and triplet projection operator as the total spin 
of the two protons in the final state, S', takes the value S' = and S' = 1, respectively. cTj 
denotes the Pauli spin matrix acting on each of the two interacting protons, i = 1 and 2. In 
terms of the Pauli spin matrices, we have Ps'=o = (1 — f?i • f?2)/4 and Ps'=i = (3 + ai • (T2)/4. 
We, then, may write, Ps'(^i ~ ^2)/2 = [((Ti — (T2) — {—)^'i(a-i x (T2)] /4, which, up to an 
irrelevant phase, is identical to the structure given in Ref. The vectors A and B in 

Eq. may be constructed from the momentum vectors available in the system, e.g., the 
relative momenta of the two protons in the initial state p and in the final state p'. Since we 
restrict ourselves to 5* and P waves for the outgoing pp system we may write 

A = ap , B = I3p' + -f{p' - 3p{p ■ p')) . (2) 

Here the amplitudes a, /3, and 7 correspond to the transitions ^Pq ^•S'qs, ^Pqs, 
and — >■ ^P2S, respectively. Note that we pulled out the linear momentum dependence, 
characteristic of P-waves, from the corresponding partial-wave amplitudes. The amplitude 
a has a strong dependence on the relative energy of the pp system in the final state refiecting 
the strong pp FSI in the state. The amplitudes (3 and 7 also depend on the relative 
energy of the pp system in the final state; however, their energy dependence is much weaker 
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than that of a due to the much weaker FSI in the ^Pq and states as compared to the 
""^^o state. 

From Eq. (jT)) exphcit expressions for any observable follow directly. E.g., we find 

^A, = .(1* X 1), , 

^^C,, = 6,,{\A\'-\B\')-2Re{A:A,) , (3) 

where Aj denotes the analyzing power and Cij the spin correlation function. Inserting the 
expressions of Eq. into Eq. (jH)) we get 

^ = \a\'+p"[\(3 + ^\' + 3x\\^\'-2Re{(3Y))] , 
—A- - 

= lap -p'^ + 7r + 3x^(171' -2Re(/57*))] , (4) 

where we introduced p'x = p' ■ p. We note that partial-wave amplitudes with even and odd 
final pp relative orbital angular momenta cannot interfere with each other due to the Pauli 
principle. 

Let us recall at this stage that the proton angular distribution seen in the experiment is 
isotropic U]. From the above equations we can see immediately that there are two possible 
scenarios for achieving such an isotropic distribution in the presence of significant pp P-wave 
contributions, namely 

1) Dominant contributions from the transitions ^Pq ^SqS and ^Sq ^PqS, but 
negligible contributions from ^P2S (7 ~ 0). 

2) Contributions from all three transitions, ^Pq — > ^Sqs, -^Sq ^PqS and ^D2 ^P2-s, 
where the latter two interfere destructively (I7P ~ 2Re(/57*)). 

Obviously, the observables given in Eq. (0} do not allow one to distinguish between the 
two scenarios and, consequently, there is no model independent way to extract the two pp P- 
wave amplitude contributions {(3 and 7) from these observables. To do that, one would need 
observables depending also on the spin of the final pp state, such as spin transfer coefficients. 
Note, however, that the combination da/dfl {C^x + 1) depends only on the amphtude a. 



Hence, a measurement of this observable would determine, in a model independent way, the 
pp S-wave contribution (a) in the final state. Similarly, a measurement of da/dQ {Cxx — 1) 
would confirm the presence of pp P-waves in the final state. It should also be stressed that 
Cxx by itself is already very interesting. Here the two angular independent terms (first two 
terms in the last equation of (@))) have opposite signs and thus tend to cancel each other. 
Consequently, this observable should be rather sensitive to the angular dependent term. 
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We now turn our attention to the results of a model for pp pprj. In Ref. 
presented a relativistic meson-exchange model for r] production in A^A^ collisions. It treats 
rj production in the Distorted Wave Born Approximation and includes both the A^A^ FSI 
and initial state interaction (ISI), the latter through the approximate procedure proposed 
in Ref. While this model yields a satisfactory description of the near-threshold cross 
section data (for pp — > pprj and for pn —>■ prirj) it fails to reproduce the recently measured 
invariant mass distributions. It should be stressed that the main objective of the present 
model calculation is not to achieve an accurate description of the existing data but rather 



to verify whether the model of Ref. 



features exhibited by the new data as discussed above. 
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In the development of a variant of the model 
cations of the vNN* vertex and the mixing parameter A in the nNN* and rjNN* vertices. 
Here, v stands for either the p- or cu-meson and A^* is the 5*11 (1535) resonance. We also use 

n 

the Paris A^A^ T-matrix [1^ as the pp FSI; the Coulomb interaction is fully accounted for as 
described in Ref. Everything else was kept unchanged. In contrast to Ref. in the 
present work we have chosen a more general gauge invariant Lagrangian jl^ for the vNN* 
coupling 



QuiNN* 



+ h.c. , 

QpNN* 



^rriN* + 
+ h.c. , 



7a. 



7/. 



niN* + 'rriN 
rriN* + fTiN 



- ida + t^oO-uiyd" 



(5a) 
(5b) 



where lj^{x), p^{x), ^ipNix) and ipN*{x) denote the u, p, nucleon and spin-1/2 A^* 
(='S'ii(1535)) resonance fields, respectively. rriN* denotes the mass of the nucleon reso- 
nance. The coupling constants Qvnn* and k„ were considered to be free parameters in the 
calculation and have been adjusted to reproduce (globally) the available data, including the 



pprj and pnr] total cross sections. The 757^ coupling at the vNN* vertex was necessary to 
achieve a reasonable fit to the data. The values obtained are: QpNN* / {^*n + = 9.5 
[/m], = -5.3 and QuNN* I ijn*N + "^a') — 6-0 [fm], = 3.8. In addition to the vNN* 
vertex given above, we have also chosen the pseudoscalar-pseudovector mixing parameter to 
be A = 0.7 at the ttNN* and rjNN* vertices jl^. All other parameter values are the same 
as given in Ref. [l^j corresponding to the case of pseudoscalar meson dominance. We refer 



to Ref. [1^ for further details of the model. 

Results for the pp invariant mass distribution based on different partial wave contributions 
are shown in Fig. [H together with the recent data at the excess energies of Q = 15 0, Q] and 
41 MeV yi . The full results are denoted by solid curves. Hereafter, these correspond to the 
calculations performed using the plane-wave basis without a partial wave decomposition and, 
as such, they include all partial waves. As is evident from the dashed curves, the observed 
peak in the region ~ (2mAr)^ is due to the strong pp FSI in the ^Sq state. The observed 
bump in the higher rripp region is largely due to the ^PqS final state (dash-dotted curves). 
The contribution from the ^P2S state (dotted curves) is very small. Contributions from 
other partial-wave states (mainly ^P2,^ -^2 -^^ Sod) are relatively small at Q = 41 MeV and 
practically negligible at Q = 15 MeV. Thus, the present model is in line with the scenario 
(1) discussed above. Overall, the shape of the pp invariant mass distribution exhibited by 
the data is nicely reproduced. However, the model tends to overestimate the data close to 
the maximum value of m^^ at Q = 41 MeV. In principle, this discrepancy might be due to 
the pr] FSI which is not explicitly accounted for in our model However, in order to 

reduce the predicted value, one needs a repulsive prj FSI which seems to be in contradiction 
with all other evidence of prj FSI effects in meson production 

BE3. 

Moreover, no such 

discrepancy is seen at Q = 15 MeV where the effect of the pt] FSI should be even larger. 
Further investigation is required to resolve this issue. 

It is important to note that the relative strength of the different partial-wave states 
depends crucially on the details of the model, and that means, specifically, on the excitation 
mechanism of the S'ii(1535) resonance in the present case. In fact, the measured pp invariant 
mass distributions can be described with the same quality as shown in Fig. ^with the ^P2S 
state contribution dominating over the ^PqS state contribution. Such a scenario can easily 
be achieved in our model by a proper adjustment of the coupling constants at the vNN* 
vertex appearing in the underlying Lagrangians (Eq. (0)). However, the resulting proton 
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angular distributions are then much more pronounced and, consequently, in disagreement 
with the experimental evidence ^. 

Although significant pp P-waves in the final state seem to be necessary for reproducing 
the pp invariant mass distribution, it is important to note that the energy dependence of 
the total cross section near the threshold energy region is basically reproduced by the pp 
FSI in the ^Sq state folded with the phase space. E.g., the model developed by V. Baru et 
al. reproduces nicely the energy dependence of the total cross section from threshold 
up to Q ~ 50 MeV with S-wave contributions alone. Results of the present model for the 
total cross section are shown in Fig. |21 Comparing the curves for the ^Sqs (long-dashed) 
and ^5*0^ +^ PqS (dash-dotted) partial waves, one realizes that the onset of the ^PqS final 
state occurs at a fairly low excess energy and that its contribution becomes increasingly 
important with increasing energy. This feature is a direct consequence of the requirement 
of reproducing the pp invariant mass distribution. However, as a result, the model now 
underpredicts significantly the data for energies close to threshold. The thin dashed curve 
corresponds to the ^Sqs contribution multiplied by an arbitrary factor of 3. This clearly 
shows that the total cross section data in the low energy region favor a larger contribution 
of the ^5*05 final state than is predicted by our model. Whether one is able to reconcile these 
seemingly contradictory properties within a consistent theoretical model remains to be seen. 
In any case, one should keep in mind that the rjN FSI, which is not included in the present 
model calculation, should enhance the ri s-wave contribution near threshold to some degree 

Q. 

Fig. ini shows the differential cross sections as a function of the rj emission angle in the 
overall cm. frame for two excess energies. The data from Ref. are basically isotropic, 
indicating a dominant 1] s-wave contribution. The theoretical results are normalized to the 
total cross section (obtained by integrating the differential cross section data) in order to 
facilitate a proper comparison of their angular dependence with the experiment. At Q = 15 
MeV the normalization factor is about 2.7, while at Q = 41 MeV, it is about 0.9. The 
dashed curves correspond to the s-wave contribution, the dash-dotted curves to the s + p 
waves, and the dotted curves to the s+p+d waves. The last are practically indistinguishable 
from the corresponding full results which are denoted by the solid lines. As expected after 
the discussion above, the angular distribution is given primarily by the s-wave contribution, 
with a small contribution from higher partial waves provided mainly by the d-wave. 
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In Fig. m the proton angular distributions in the overall cm. frame are shown together 
with the data from Ref. j6|. Here again the model predictions are normalized to the data 
using the same factor mentioned above. Evidently, at Q = 15 MeV the result is isotropic; it 
is dominated by the ^SqS (dashed curve) state followed by the ^PqS (dash-dotted curve) state. 
Contributions from other partial waves are practically negligible. At Q = 41 MeV there is 
a small contribution from the ^P2S (dotted curve). Obviously, the destructive interference 
with the ^Pqs state canceling the angular dependence (see Eqs.(jH)) is incomplete resulting 
in a noticeable angular dependence which, however, is still compatible with the experiment. 
The difference between the dotted and solid curves is due to a small contribution from the 
— >^ D2S plus ^P2,^ F2 — >^ SqcL amplitudes. 

Fig. El shows the prediction for the analyzing power as a function of the r] emission angle 
in the overall cm. frame. The data are from Ref. P]. The dashed curves correspond to 
the 7] s + p wave contributions while the dash-dotted curves to the s + p + d waves. The 
solid lines are the full results. We note that the s-wave contribution alone yields Ay = 
(see Eqs.(@])), so that any non- vanishing result must necessarily involve higher partial waves. 
Furthermore, judging from the shape exhibited by the analyzing power, the present model 
yields a vector meson dominance over the pseudoscalar meson in the excitation mechanism 
of the 5*11(1535) resonance as discussed in Ref. Although the data indicate some 

contribution from partial waves higher than the s-wave, they are not sufficiently accurate 
to make a definitive statement as to the size of their contribution. 

In Fig. ini we present predictions for the spin correlation function C^x at Q = 41 MeV 
as a function of the proton angle in the overall cm. frame. As can be seen from Eq. (^, 
the ^Pq ^SqS partial wave alone leads to a constant value of C^x = 1- Adding the 
^5*0 —>■ ^Pqs contribution (dashed curve) yields a small value of Cxx ~ —0.3. Including 
also the ^D2 ^P2S contribution one obtains the result represented by the dotted fine. 
(In this context note that the ^Sq — * ^Pqs as well as the ^D2 — > ^P2S contributions alone 
would give rise to Cxx = —1, cf. Eq. (jlj)). Other partial- wave contributions do not change 
Cxx qualitatively as is evident from the full result (solid curve). Thus, in our model, the 
cancellation of the \a\^ term and the 1(3 + 7P in Eqs. Q is almost complete! This strongly 
enhances the relative importance of the angular dependent term in Cxx- Note that Cxx = Cyy 
for all of the three partial- wave contributions discussed explicitly above (cf. Eqs. 0). 

In Fig. [3 predictions for the pr] invariant mass distribution are shown together with the 
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data '6]. Again, the dominant contributions are from the ^5*05 (dashed curves) and ^Pqs 
(dot-dashed curves) final states; the contribution from the ^P2S state (dotted curves) is 
neghgible. At Q = 41 MeV one sees also some effects from other partial-wave states arising 
mainly from the '^^2;'^ ^2 SqcI amplitude. The overall shape of the measured invariant 
mass distribution is reproduced. The observed discrepancies in the details, especially at 
Q = 41 MeV, are not easy to understand in view of the nice agreement between calculated 
and measured pp invariant mass distributions. 

Summarizing our results, we have shown that the currently available data on rj production 
in pp collisions near the threshold energy can be largely understood in terms of a few S- 
and P-wave amplitudes. For a completely model-independent extraction of the relevant 
amplitudes, however, observables independent of those presently available are required. In 
this connection, the spin correlation function, either Cyy or C^x, is suited to further constrain 
the ^PqS and ^P2S final state contributions. In any case, the final pp P-wave contribution 
is crucial for explaining the measured invariant pp mass distribution, especially, aX Q = 
AlMeV. Our model calculations show that, the dominant amplitudes are ^Pq ^Sqs and 
^5*0 — * ^PqS- It should be stressed that in order to quantify the role of the rjN interaction 
in pp — ^ pprj it is important to first understand the role of higher A^A^ partial waves. 

Finally, we note that the present work illustrates the possibility of using meson production 
processes in A^A^ collisions to study the properties of nucleon resonances in terms of a few 
partial-wave amplitudes. In particular, the present model prediction for the relevant partial- 
wave amplitudes depends very sensitively on the details of the model and especially to the 
excitation mechanism of the 5*11(1535) resonance. This offers an excellent opportunity to 
study some of the properties of the 5*11(1535) resonance using the r] meson production 
reaction in A^A^ collisions which would not be possible to investigate in more basic reactions 
such as7 + A^^?7 + A^ and M + N ^ r] + N. 
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FIG. 1: Final pp invariant mass distribution in pp pprj as a function of invariant mass squared, 
nipp, at an excess energy of Q = 15 MeV (upper panel) and Q = 41 MeV (lower panel). The 
dashed (dash-dotted) curves correspond to the ^Sqs (}Sos +^ Pqs) final state contribution. The 
dotted curves correspond to ^5*05 -|-^ Pqs +^ P2S; its is indistinguishable from the solid curve in the 



upper panel. The solid curves are the full results. The data are from Ref. 



(filled circle) and 



3rps 

Ref. [3] (open circle). The latter have been normalized by an arbitrary factor of 0.66 in order to 
facilitate the comparison of the shape with the former data and the present model prediction. 
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FIG. 2: Total cross section for the reaction pp — > pprj as a function of the excess energy Q. The 
solid curves represent the full results. The dashed (dash-dotted) curve corresponds to the ^Sqs 
{^Sqs +^ Pqs) final-state contribution and the dotted curve to the ^Sqs +^ Pqs +^ P2S contribution. 
The thin dashed curve is the ^Sqs contribution multiplied by an arbitrary factor of 3. The data 
are from Ref. y|. 
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FIG. 3: Angular distribution of the emitted rj meson in the cm. frame of the total system at 
an excess energy of Q = 15 MeV (upper panel) and Q = 41 MeV (lower panel). The dashed 
(dash-dotted) curves correspond to the r] meson s-wave {s + p-wave) contribution. The dotted 
curves correspond to the s + p + d-wave contributions which are practically indistinguishable from 
the corresponding full results represented by solid curves. The data are from Refs. P, ^• 
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FIG. 4: Angular distribution of the proton in the final state in the cm. frame of the total system 
at an excess energy of Q = 15 MeV (upper panel) and Q = 41 MeV (lower panel). The dashed 
(dash-dotted) curves correspond to the ^Sqs (^Sqs Pqs) final-state contribution. The dotted 
curve corresponds to the ^Sqs +^ Pqs P2S contribution and the solid curve represents the full 
result. The data are from Ref. ;G]. 
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FIG. 5: Analyzing power for the reaction pp — > ppr] as a function of rj emission angle in the cm. 
frame of the total system at an excess energy of Q = 15 MeV (upper panel) and Q = 41 MeV 
(lower panel). The dashed (dash-dotted) curves correspond to the partial waves contributions with 
/' < 1 (/' < 2). The solid curves represent the full results. The data are from Ref. 
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FIG. 6: Predicted spin correlation function C^x for the reaction -pp pprj as a function of final 
proton angle in the overall cm. frame at an excess energy of Q = 41 MeV. The dash (dash-dotted) 

curve corresponds to the ^Sqs Sqs +^ Pqs) final state contribution. The dotted curve corresponds 
to the ^Sqs +^ Pqs +^ P2S contribution and the solid curve represents the full result. 
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FIG. 7: Same as Fig. ^ but for the pr] invariant mass distribution. The data from Ref. (open 
circle) have been normaUzed by an arbitrary factor of 0.66 in order to facihtate the comparison of 
the shape with the data from Ref. 6] and the present model prediction. 
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